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I. INTRODUCTION 


In the last years much effort was put, on both experimental and theoretical sides, into 
understanding the spin structure of the proton (for reviews of the EMC spin crisis, see |1J). 
New data for the polarized deep inelastic scattering (DIS) structure function g 1 from SLAC 
[2J and the SMC [d] are consistent with the previous EMC J3J result. Two different scenarios 
can explain the data, large anomalous sea-quark or large anomalous gluon interpretation 
M which are connected with the problem of the anomalous gluon contribution to (j \. In 
fact, a renormalization group transformation can shift the anomalous contribution to the 
quark sector or to the gluon sector. In this way, either the anomaly manifested in the matrix 
element of the axial quark current is responsible for the large polarization of the sea or it can 
be perturbatively taken into account in the gluon coefficient function, while the quark matrix 
element is anomaly-free. The dependence of this renormalization group transformation on 
the gluon momentum fraction x is still controversial. 

From the physical point of view, it is more natural that the anomaly is attributed to the 
gluons. The gluonic coefficient function then absorbs only the short-distance contributions, 
while the quark distribution corresponds to a conserved operator and may be invoked in 
the low-energy description of the proton. The corresponding picture of the nucleon spin 
structure is most popular and the gluon polarization is considered to be its most important 
unknown ingredient. The problem of its x-dependence is the “physical” counterpart of the 
above-mentioned renormalization group transformation. 

The accounting for the anomaly by the finite renormalization transformation is effectively 
resulting in the substitution in the partonic expression for the structure function g \, 

Qi 

Aq ^ Aq = Aq- ^-Ag, (1.1) 

for the first moment of the spin-dependent quark distribution for each flavour. The over¬ 
simplified “naive” x-dependent analogue is just 

ry 

A q(x) = A q(x) — -^AG(x). (1.2) 
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However, such an expression is not compatible, in principle, with the fact that the gluon 
contribution starts at one-loop level only. The simplest, consistent approach was explored 
in Ref. ]8| a few years ago, by taking the infrared (IR) finite part Sir, of the box diagram, 
responsible for the photon-gluon interaction: 

Aq(x) = A q(x) - ^ J* yS IR Q A g(z). (1.3) 

However, the choice of the hnite part is ambiguous and the specific role of the axial 
anomaly remains unclear. This is why the attempt followed in Ref. |J to perform the 
decomposition of the box diagram, referring to its kinematical structure. The anomaly then 
contributes to the singlet structure, associated with the spin structure function gr = 9i + 92 
of the photon-gluon scattering^]: 


9t = 


-A g n . 


(1.4) 


2 n{n + 1) 

Attributing the anomalous gluon contribution to this structure function resulted in its 
1 — ^-dependence. Later, this expression was obtained and exploited in different approaches 
mill ■ However, both derivations and interpretations of this result do not seem complete 


and some more solid ground is desirable. 

The anomalous gluon contribution plays an exceptional role in the general classification of 
perturbative QCD contributions. Formally, it is a part of the next-to-leading order (NLO) 
contribution. However, due to the well-known growth of the first moment of the spin- 
dependent gluon distribution, compensating one power of a s makes it essential also at leading 
order. The recently calculated two-loop, spin-dependent anomalous dimensions [0)01 niake 


it now possible to perform the complete NLO analysis of the experimental data [15-18 


These NLO calculations were performed in the dimensional regularization using the ’t Hooft- 
Veltman-Breitenlohner-Maison (HVBM) scheme [O] in which y 5 does not anticommute in 


1 The recent analysis [jlCj] confirms that the zero first moment of this structure function is mani¬ 
fested only for regularization schemes, providing the zero moment of g\ as well. 


3 























the unphysical space-time dimensions. In this minimal subtraction (MS) scheme both chiral 
invariance in the non-singlet sector and the one-loop character of the singlet axial anomaly 
are explicitly broken and their restoration requires an additional finite renormalization. 


As the HVBM scheme is used, a further renormalization group transformation is per¬ 
formed [T5LB0| in order to make the result compatible with the standard factorization pre¬ 
scription and the low-energy intuitive description of the proton. However, only the value of 
the first moment was fixed in that procedure, so that there remains an ambiguity in this 
transformation. Note that the function 1 — x was also proposed for this purpose (jl!| . 

In the present article, we are suggesting the non-local generalization of the axial anomaly 
based on the canonical Ward identities for light-ray operators. This allows us to give the 
natural description of both anomaly-free singlet quark distribution and anomalous gluon 
contribution, which fix the rc-dependence of the renormalization group transformation. The 
final result is a rigorous proof of the 1 — x-behaviour of the anomalous contribution in NLO. 
Here, we do not address the issue of higher loop corrections to the generalized axial anomaly. 

The paper is organized as follows. Section II discusses the chiral invariance breaking 
of flavour non-singlet light-ray operators due to the renormalization, and their restoration 
by a finite renormalization. Then we derive Ward identities for light-ray operators and 
compute the non-local singlet axial anomaly, which can be expressed as a divergence from a 
generalized topological current. In Section III we use our results to define the anomaly-free 
singlet quark distribution and show the consequences for coefficient functions and evolution 
kernels in NLO approximation. 


II. AXIAL ANOMALY OF LIGHT-RAY OPERATORS 

As mentioned before, chiral invariance is broken in the HVBM scheme; however, as it is 
known from the renormalization of the axial current, it can be restored by a finite renormal¬ 
ization, so that the non-singlet Ward identity will be fulfilled. This finite renormalization 
constant z NS can be computed for massless QCD from the requirement that the anticom- 
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mutativity of y 5 is effectively restored |2l 


, i.e. 


„NS 


■5,a 
Jfi 


W>) = 


CM. a 


with z Nb = 1- —C F + 0(a 2 s ), 

' 77 


( 2 . 1 ) 


where j^ ,a = ^[y^,y 5 ]A a, 0 with y 5 = hy^^y 3 and j“ = T/jy^A 0 ^ are the axial vector 
and the vector current, respectively, A a is a flavour matrix, the symbol [...] denotes minimal 
subtraction, and C F — | is the usual QCD colour factor. Flavour and colour indices are 
suppressed for simplicity. 

The one-loop character of the singlet axial anomaly can also be ensured by some finite 
renormalization of the current. This required a two-loop calculation of the Ward identities 
sandwiched between the two-gluon state. The a s correction to the corresponding z factor, 
z s = 1 — 9 ^C F + 0(a 2 ), coincides with the non-singlet factor z NS . Because of the appearance 
of so-called light-to-light subdiagrams, this coincidence is spoiled beyond the one-loop level. 
Note that the one-loop approximation of z s can also be calculated by the requirement that 
the singlet Ward identities are fulfilled for the two-quark state. Since at leading order the 
axial anomaly does not appear in this Ward identity, it follows that the finite renormalization 
factor is the same as that for the non-singlet channel. 

The same problems as discussed above also occur for composite operators on the light- 
cone, which appear in the definition of spin-dependent quark distribution functions. For 
technical reasons we start with a more general definition of bi-local operators, which are not 
necessarily on the light-cone: 


0 5 ' a {x,y) = ^(x)U 8 (x,y) 




A y(y), 


( 2 . 2 ) 


U a (x,y ) = P exp | —ig drA^xr + y[l - r})t k {x^ - j/^j . 


Here, U s (x,y) ensures gauge invariance, where the gauge held A% is path-ordered along 
a straight line connecting the fermion fields. For light-like distances, i.e. (x — y) 2 = 0, 
these operators can be expanded in terms of local twist-2 and twist-3 operators. We set 
x = Kix,y = k 2 x, where x is a light-cone vector and, after contraction with we get the 


5 











leading twist-2 light-ray operators ||22||p : 


0 5 ’ a (n i, K 2 ] X) = X^O^ a (KiX, k 2 x ), 

= ^(k,ix)U(kiX, k 2 x) ftrf\ a '4’{K 2 x). (2.3) 


To determine the finite non-singlet renormalization constant in the forward case, we 
require (for massless QCD), in analogy to Eq. flO), the validity of 


dx z Nb (x) ( 0 5,a ( 0,kx;x) i/iip) = ([O a ( 0, k; x)]-0-0) y 5 , 


(2.4) 


where O a is analogous to the definition of 0 5,a in Eq. (|2.3|) , but without y 5 matrix. I 11 the 
dimensional regularization using the HVBM scheme, the leading-order result (restricted to 
the forward case) is 


z NS (x) = 5(1 — x) — —2C(f(1 — x) + O(ay), (2.5) 

7r 

while in the Pauli-Villars regularization, chiral invariance holds true without finite renor¬ 
malization. 

As in the case of the axial current one expects, in leading order, that the same finite 
renormalization as in Eq. ( p.5|) has to be performed for the singlet light-ray operator. Beyond 
the leading order this is no longer true and the question is: How can we fix this finite 
renormalization constant? 

After we have seen that the restoration of chiral invariance requires an x-dependent 
finite renormalization, a second question arises: Does the axial anomaly also depend on x? 
To answer this question we use the equation of motion to derive Ward identities for the 
divergence of the non-local operators ( [2.2| ). A straightforward calculation provides that the 
divergence of this operator is 


2 Since x is an external four-vector that can be kept four dimensional in the dimensional- 
regularized operator vertex, it follows that in the HVBM scheme the relation [^,y 5 ] = 2 ^y 5 


is valid, so that indeed Eq. (T3) comes from the definition (2.2). 
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(<9x + %) Ol’ a (x, y) = fi EOM (x, y) + i^(x)U s {x, y)^m a ^{y) - 

ig f dr ^(x)H s (x, xr)y a F Q/3 (xr + y[ 1 - r])(xy - yy) x 

J 0 

U s (y[ 1 - t], y)^ 5 \y(y), ( 2 . 6 ) 

fi EOM (x, ?/) = $(x)([p(x) - im]U s (x,y) + U s (x,y)[p(y) - im^X^y), (2.7) 

where fi EOM (x, y) denotes the equation of motion operators, F a y = F£pt a is the held strength 
tensor, and m“- = (ny + m,j) A“- is a mass matrix. 

For (x — y ) 2 7 ^ 0 the operator is ultraviolet-finite and the Ward identity ( ^ 6 |) should be 
satisfied. However, the situation will be changed if we go to a light-like (not only short, as 
usually discussed) distance. Then the operator has to be renormalized, and anomalous terms 
can occur. For instance, if we naively anticommute P(y) with y 5 in the second expression 
in the r.li.s. of Eq. ( |2.7| ) then the use of the equation of motion will provide only contact 
terms. However, in the HVBM scheme the non-anticommutativity of y 5 gives an anomalous 
term, which should be cancelled by a finite renormalization (see the discussion in [[hj]), so 
that the anticommutativity of y 5 is effectively restored as discussed above. We calculated 
this anomalous term and after an appropriate definition of the appearing operators the finite 
renormalization constant (|2.5|) was extracted from the Ward identity. As discussed above 
for the local case, because of the absence of the singlet axial anomaly in the Ward identities 
for the Green functions of quark fields at one-loop order, it follows that z s (x) coincides with 
z NS (x) at this order. 

For the singlet case (a = 0, = 5 t]1 rri^ = 2 inAj), we expect an anomalous term of the 

form e a/3lS F a g(- ■ -)F 7 s(- • •), which can be computed from the difference of the l.h.s. and r.h.s. 
of the Ward identity (|2.6| ) sandwiched between the two gluon states (see diagrams in Fig. 
1). We now set x = x,y = K2X with x 2 = 0 and choose the light-cone gauge, i.e. xA = 0. 
After a straightforward calculation, it turns out that in both Pauli-Villars regularization 
and HVBM scheme the same anomaly appears, so that the singlet Ward identity is actually 
given by 
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/ \ Oi /■! rl—x 

(%!*) + ^ 2 X)) ^) = « 2 ; 5) + ^iV/ J o dx y o d|/2x 

F lv (M 1 - x) + k 2 x\x) F a ^ u ([k 2 ( 1 -y) + «i y]x) , (2.8) 

where 

Of(ki, k 2 ;x) — —ig f dr '^(Kix) r j a F al3 (rx)xp n / 5 'ilj(K 2 x). (2.9) 

J H\ 

Here, F a ^ = \e^F^ with e 0 i 2 3 = +1, a s = fy, and Nf is the number of flavours. Here 
we applied the equation of motion and neglected the quark masses. In the local case, which 
follows from setting K\ = k 2 = k, the operator 0^°(kx, kx) coincides with the local singlet 
axial current j^ Q (nx) and from Eq. ( |2.8|) we recover the well-known expression for the local 
anomaly: j^NfF“ u (kx) F afJ,u (kx). 

Finally, we introduce a non-local generalization of the topological current 


(X.C 


K»(x,y) = j^N I e^ 


K( x ) d i3 A °(y) - y^K( x ) Ab D(v) A y(v) + l x <-> y} 


( 2 . 10 ) 


with x = K\X and y = k 2 x, so that the anomaly in Eq. (|2.8|) can be written as a divergence 


(9(«,j) + 9f, 2 s>) K f (n ii, n 2 i) = ^ N f F( k 2 x) 


( 2 . 11 ) 


We are now able to dehne the anomaly-free operator 


O 5 ^ 0 (kix, k 2 x) = 0^°(kiX, k 2 x) - 

2 [ dx [ dy K tl ([K 1 (l^x) + K 2 x]x,[K 2 (l~y) + K 1 y]x). (2.12) 

Jo Jo 


III. CONSEQUENCES FOR POLARIZED DISTRIBUTION FUNCTIONS 

As shown in the previous Section light-ray operators possess anomalous contributions. 
For the non-singlet case these anomalies are a pure artefact of choosing a non-invariant 
chiral renormalization scheme. The polarized quark distribution functions should be defined 
in a chiral invariant manner, so that in a general renormalization scheme an additional finite 
renormalization is necessary 
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Aq NS (x,Q 2 ) = £j$ s (v) I ^^y(-P,S|[O ! " r ' s (0, K ;i)]|P,S)e i < 1:/ »>“ (ip >, (3,1) 

where S p denotes the polarization vector of the nucleon and the renormalization point square 
y 2 is set equal to the momentum transfer squared Q 2 . The leading-order approximation of 
z 5 (y) HVBM scheme is given in Eq. (|2.5| ). This hnite renormalization affects the NLO 

approximation of both quark coefficient functions and splitting kernels and agrees with the 


additional renormalization group transformation performed in the NLO calculation []T3|,[T4 
Because of the anomaly, the naive definition of the singlet distribution function 


A E(x,Q 2 ) = Aqi(x,Q 2 ) + Aqi(x,Q 2 ), Aq^x.Q 2 ) = Aq^-x.Q 2 ) 

i=u,d 

dK (p,S|[O 5 ’°(0,k;x)] +{k -»■ -K}\P,S)e i{ - x/y )^ p ) (3.2) 


/ 


2i t(xS) 


(please note that AS is defined in terms of the operator O 5,0 ) cannot be interpreted as 
probability for Ending a polarized quark flavour singlet configuration with given longitudinal 
momentum fraction x. In addition to the hnite multiplicative renormalization^], it is also 
necessary to remove the axial anomaly from the definition of this function [5-]7[]. Thus, one 
has to define the singlet distribution function in terms of the anomaly-free operator ( p,12|) , 
which provides 


AE(x, Q 2 ) = AE(x, Q 2 ) - k(x, Q 2 ), (3.3) 

where 

k(x,Q 2 ) = 2j^(l-y) J (PS|i'‘P,(0, K i) + { K ^- K }|pS)e i (* / <'l'‘ ( * p) , (3.4) 

and the generalized topological current K M is dehned in Eq. ( p,10|) . 

The problem that K^ is gauge-variant, and thus that k(x, Q 2 ) contains also unphysical 
components, can somehow be resolved by the choice of a physical gauge. In the light-cone 


3 This renormalization is due to the diagrams without two-gluon intermediate states, which are 
the same for the non-singlet and the singlet case, so that z^(y) = Zjh s (y), as discussed in Section 

II. 
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gauge, the gauge-invariant twist-2 gluon operator 


G(k!,k 2 ;x) = ix a F aa9 (K 1 x)F^(K 2 x)x^ 


(3.5) 


can be expressed in the forward case by 

G( 0, k;x) = - ixd^K^O, 


KX). 


(3.6) 


Furthermore, from the definition of the gluon distribution function 


Ag(x,Q 2 ) = 


dn 


(. P , S |G(0, k; x) + {k -4| P, S') 4 K(iP) , (3.7) 


x(xP) J 2n(xS) 

and from Eqs. (|3.4|) and ( |3.6| ), and after performing a partial integration, we find that 
k(x, Q 2 ) is actually given in terms of the gauge-invariant gluon distribution function. Thus, 
to remove the axial anomaly of the polarized quark distribution function it is sufficient to 
subtract a certain amount of the gluon distribution function: 

AE(x, Q 2 ) = AE(x, Q 2 ) — K (x) 0 Ag(x, Q 2 ), K{x) — - -Nr(l—x), (3.8) 

7r 


where the convolution is dehned as 


A(x) <S> B(x) — [ dy [ dz 5{x — yz)A{y)B{z). 
Jo Jo 


(3.9) 


Indeed, removing the axial anomaly is equivalent to the following (additive) renormal¬ 
ization group transformation: 


A C g (x) = A C g (x) + K ® AC q (x), AC g (x) = AC g (x), (3.10) 

^Pgg ( x ) = ^Pgg( x ) + di ® A P gg (x), AP gq {x) = AP gq (x), 

AP qq (x) = AP qq (x) - K ® AP gq (x), 

A Pqg(x) = A P qg (x) - ~ K ( X ) + K ® [APgg - APgg ~ A Pgg <g) K](x), 

where A4 are the coefficient functions, AP ig are the spin-dependent splitting kernels, and 
/3 = fi-jj- L g(g) is the renormalization group coefficient of the running coupling constant. 
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Because of gauge invariance the operator product analysis suggests that the zero moment 
of A C g vanishes. After the transformation (p. 1C|) is performed, the gauge-variant axial 


anomaly contributes to the gluonic sector, so that the zero moment is now given by 

J dx A C g (x) = J dx K{x ) + O (a^j = Nf + O (af) . (3-11) 

For completeness we give the difference of the splitting kernels in NLO: 

A P gg {x) - AP gg (x) = - (A P qq (x) - A P qq (x)) 

= (^tt) 2( -' f -^ f ~ x ) + ( 2 + x ) In x ], 

A Pgg(x) - A P qg (x) = N f {C F [(1 - x)(l - 41n(l - x) + 21nx)] + 

Ca [(1 — x)(—16 + 41n(l — x)) — 4(2 + x) lnx] }, (3.12) 

where the Casimir operator Ca is equal to the number of colours. Note, that these differences 
will vanishes in the limit x —> 1. For small x they contribute only to the subleading 
behaviour of the NLO result (the leading terms are given by ln 2 x). 

Another opportunity to perform the evolution coincides with the one proposed by Cheng 
. Namely, one should perform the evolution in the gauge-invariant (say, MS) scheme and 
afterwards restore the anomaly-free distribution by applying Eq. (|3.8|) . The 1 — x behaviour 
in this approach is actually coming from the mass term in the box graph. The contact with 
our derivation may be achieved by noting that the cancellation of normal and anomalous 
divergence, resulting in the effective conservation of axial current in the limit of infinite 
quark mass, is valid for the non-local anomaly as well. It is especially clear for the Pauli- 
Villars regularization, when the contribution of regulator fermions (calculated in Section II) 
is looking, up to the sign, precisely such us that of the quark masses, which is the starting 
point of the approach of Cheng. 

For practical purposes, irrespective of the used evolution scheme, it is possible to define 
an effective gluon distribution, which is just the combination appearing in rq, i.e. the 
convolution of (1 — x) with A g(x): 

A g cS (x, Q 2 ) = 2(1 — x) ® A g(x, Q 2 ), (3.13) 
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so that the first moments of effective gluon distribution coincide with the “original” one, 
while the structure function gq has at leading order the simple partonic form, suitable for 
the extraction of partonic distributions from the experimental data |24|] : 


gi{x,Q 2 ) = ^A t(x,Q 2 ) = i 


A E(;r,Q 2 


Os 

2n 


N f Ag ett (x,Q 2 ) 


(3.14) 


It is convenient to have the evolution equation directly for effective distribution functions. 
While the diagonal kernels are not changed, the moments of the off-diagonal kernels are 
changed in a straightforward manner: 

2 


, „ AP w (n), 


(«) = A±h APm(n) , 


(3.15) 


n(n+ 1 ) —19--' 2 

Note, especially, that the influence of the effective gluon distribution to the quark evolution, 
governed by A P^, appears to be much less singular in n. We are not presenting the explicit 
form for effective NLO anomalous dimensions, which are rather lengthy. At leading order, 


when the x dependence is easily restored, Eq. (|3.15|) results in the equations 
AP s f (x) = ^C f [3(x - 1 ) - (x + 2) In f], A Pg = i - x). 


(3.16) 


The first moments of the splitting kernels provide an important check for the normal¬ 
ization of the non-local contribution and determine the evolution in the corresponding sum- 
rules. So we summarize the consequences for the zero moment of the splitting kernels, which 
come from general renormalization arguments of the axial vector current and the topological 
current, which are verified up to two- and three-loop order, respectively |^5||. From current 
conservation of j^° = — k M (here, j^’° refers to the original definition in terms of quark 

fields), the Adler-Bardeen theorem, and the renormalization properties of gauge-invariant 
operators, it follows that 

[ dx A Pgg{x) = f dx A Pgg(x) + "fj = —- + , 

J 0 J 0 ^ 


dx A Pgq(x) = 


dx A Pgq(x) = 


7 .i 


JV/[a./( 2 tt)]’ 


/ 0 


dx A P qg (x) = dx A P qg (x) = 0 , 
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[ dx A P gq (x) = f dx A P gq (x) - 7 j = 0 , 
/ 0 '-'0 


(3.17) 
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where the three-loop order approximation for the anomalous dimension of the axial vector 
current 7 j and for the /3-function are knownQ |25|j26j : 


7 j = 


a l\ 

Att) 

i i—C A 

g 47T V 3 


(3.18) 


6C F N f + (^) [~C A + + 18C f ) C f N,, 

§*')-(£)’( T C, *- 2C > N '- 1 T C * N >)- 

(s ) 3 (lr cl + cpN > - ii c ^ N f ~ + i c * N f) ■ < 3 - 19 > 


These results allow [^7] the extraction of the anomaly equation renormalization for a 
number of loops exceeding that of 7 j by 1. For the leading three-loop contribution, this 
coincides with the calculation of Anselm and Johansen p 8 |, while the four-loop correction 
require a one-loop finite correction to the gluon matrix element of the topological current 
29|. As a result, the four-loop correction takes the form 

3 fn \ 4 / 214 4 


Zi = - 


(s) 6CpN > + (s) + + 18Cf ) CfN >■ (3 - 20) 
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Figure 1: Diagrammatical representation for tlie l.h.s. (a) and for the r.h.s. 
(b) of the non-local singlet Ward identity (2.8) in the one-loop approximation. 
The symbol dO refers to the divergence + d^ K2 x)) K 2 %) and 

Of denotes the operator Of(Ki,K 2 ;i’) clehned in Eep (2.9). 
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